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NOTE ON THE PICARD METHOD OF SUCCESSIVE APPROXIMATIONS. 

By Dunham Jackson. 

The Picard method of successive approximations, as applied to the 
proof of the existence of a solution of a differential equation of the first 
order, is commonly introduced somewhat after the following manner: 

" We shall develop the method on an equation of the first order 

(1) !-/(*,»). 

supposing first that the variables are real. We shall assume that the 
function / is continuous when x varies from x to x + a and when y 
varies between the limits (y — b, j/ + b) ; that the absolute value of the 
function / remains less than a positive number M when the variables x, y 
remain within the preceding limits ; and, finally, that there exists a positive 
number A such that we have 

\f(x, V) -f(x,y')\< A\y -y'\ 

for any positions of the points (x, y) and (x, y') in the preceding region. 

" Let us suppose, for ease in the reasoning, a > 0, and let h be the 
smaller of the two positive numbers a, b/M. We shall prove that the 
equation (1) has an integral which is continuous in the interval (x , x + h) 
and which takes on the value y for x = x ." 

This particular language is quoted substantially from Goursat's 
Mathematical Analysis, translated by Hedrick and Dunkel,* except 
that the statement there is for a pair of differential equations in two un- 
known functions. The italics are kept from the original French. After 
the proof has been given, the following remark is added :f 

" If ... we go over the proof again, we see that the condition h < b/M 
is needed only to make sure that the intermediate functions y u y«, ■ ■ • 
[[the successive approximations to the solution] do not get out of the 
interval (y a — b,y + 6), so that the functions /(x, yi) shall be continuous 
functions of x between x and x + h. If the function /(x, y) remains 
continuous when x varies from x to x + a, and when y varies from 
— so to + t», it is unnecessary to make this requirement." 

* Vol. 2, part 2, pp. 61-62. 

t Loc. oit., p. 64; the statement is again simplified from two unknowns to one, in quoting. 
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76 DUNHAM JACKSON. 

The purpose of this note is to point out that even if f(x, y) is originally 
defined only in a rectangle, it is a simple matter to extend its definition 
outside the rectangle, so that the conditions of the hypothesis shall hold 
for x si x S Xo + a and for all real values of y, the Lipschitz condition 
as well as the mere continuity. It is sufficient, for example, to let 

fix, y) = f{x, y + b), y S y + b, 
fix, y) = fix, 1/0-&), y ^ y - b. 

The process of successive approximations then gives, at a single stroke, 
a function y(x) which is defined and satisfies the differential equation, 
with the extended definition of fix, y), for x S x ^ x + a. It satisfies 
the original equation as long as the x and y of the solution remain within 
the original rectangle, whatever the behavior of the approximating 
functions may be. The solution is unique as long as it stays in the 
rectangle. The original equation of course has no authority outside its 
own domain, and corresponding to the infinitely many possible ways of 
extending the definition of / there will be infinitely many different ex- 
tensions of the solution, if it leaves the rectangle before x reaches x 4- a. 

It may seem that this observation is trivial, and it is perhaps hardly 
probable that it is made here for the first time;* but its" omission from 
standard presentations of the subject is notable. In the treatise ahead y 
quoted, for example, after the Cauchy-Lipschitz proof has been explained, 
the two demonstrations are compared as follows :f 

" Cauchy's first method [the Cauchy-Lipschitz method] and that of 
the successive approximations give, as we see, the same limit for the 
interval in which the integral surely exists. But from a theoretical point 
of view Cauchy's method is unquestionably superior: we shall show, in 
fact, that this method enables us to find the integral in every finite interval 
in which the integral is continuous." 

Again, the Encyclopedic des sciences mathematiques, in the article 
Existence de Vintegrale generate. Determination d'une integrate particuliere 
par ses valeurs initiates, after going into some detail on the question of the 
length of the interval of convergence, says: J 

" On ne connatt encore aucun moyen de determiner l'intervalle exact 
dans lequel la methode de E. Picard converge. Suivant les cas, cet 
intervalle peut embrasser, comme dans la methode de Cauchy-Lipschitz, 
tout l'intervalle de r^gularite de la solution, ou etre au contraire plus 
petit que l'intervalle de convergence des series de Taylor en (x — x ) qui 
repr^sentent la solution, quand elle est holomorphe." 

* Since this note was written, I have learned that Professor Wedderburn made essentially the 
same suggestion, in unpublished form, a number of years ago. 
t Loc. cit., p. 73. 
% Tome 2, volume 3, fascicule 1, pp. 14-15. 
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While the present remarks do not perhaps invalidate either of these 
statements, it does seem fair to say that they have a bearing on the 
comparison. 

It is readily seen that the method outlined above can be extended to 
the case of a system of n differential equations in n unknown functions. 
If there are two equations, for example, 

dy dz . . 

Jk = K x > y> 2 )> dx = *( x ' y ' z >' 

with right-hand members defined for x = x ^ x + a, Vo — b ^ y 
= 2/o + b, zo — c s z S zo + c, it is possible to set 

f(x, y, z) = fix, yo + b, z), y S y + b, z — c S z ^ z + c; 

fix, y, z) = fix, y, zo + c), y - b ^ y ^ y + b, z S Zo + c; 

fix, y, z) = fix, 2/0 + b, z Q + c), y^yo + b, z s z + c; 

and similarly in the other regions of the yz-plane, with a corresponding 

treatment for <p. More concisely, for any number of dimensions, the 

value of each function at any point outside its original domain of definition 

is to be the same as the value which it has at the nearest point of that 

domain. 

The same method, though of course not the same formulas, can be 

used even if the original domain is not rectangular, provided that it has 

a moderately regular boundary, so that the functions can be extended 

across the boundary with the requisite degree of continuity. 

The University op Minnesota, 
Minneapolis, Minn. 



